Transport of interacting electrons through a double barrier in quantum wires 
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We generalize the fermionic renormalization group method to describe analytically transport 
through a double barrier structure in a one-dimensional system. Focusing on the case of weakly 
interacting electrons, we investigate thoroughly the dependence of the conductance on the strength 
and the shape of the double barrier for arbitrary temperature T. Our approach allows us to system- 
atically analyze the contributions to renormalized scattering amplitudes from different characteristic 
scales absent in the case of a single impurity, without restricting the consideration to the model of 
a single resonant level. Both a sequential resonant tunneling for high T and a resonant transmis- 
sion for T smaller than the resonance width are studied within the unified treatment of transport 
through strong barriers. For weak barriers, we show that two different regimes are possible. Mod- 
erately weak impurities may get strong due to a renormalization by interacting electrons, so that 
transport is described in terms of theory for initially strong barriers. The renormalization of very 
weak impurities does not yield any peak in the transmission probability; however, remarkably, the 
interaction gives rise to a sharp peak in the conductance, provided asymmetry is not too high. 

PACS numbers: 71.10.Pm, 73.21.-b, 73.23.Hk, 73.63.-b 



I. INTRODUCTION 

Effects related to the Coulomb interaction between 
electrons become increasingly prominent in systems of 
lower spatial dimensionality as their size is made smaller. 
Recent experimental progress in controlled preparation of 
nanoscale devices has led to a revival of interest in the 
transport properties of one-dimensional (ID) quantum 
wires. Owing to the particular geometry of the Fermi 
surface, systems of dimensionality one are unique in that 
the Coulomb correlations in ID change a noninteracting 
picture completely and thus play a pivotal role in low- 
temperature physics. 

A remarkable example of a correlated ID electron 
phase is the Luttinger-liquid model (for a review see, 
e.g., Refs. EBI)- In this model, arbitrarily weak interac- 
tions ruin the conventional Fermi liquid phenomenology 
by essentially modifying low-energy excitations across the 
Fermi surface. As a result, the tunneling density of states 
develops power-law singularities on the Fermi surface. 
Moreover, interactions between oppositely moving elec- 
trons generate charge- and spin-density wave correlations 
that lead to striking transport properties of a Luttinger 
liquid in the presence of impurities. In particular, even a 
single impurity yields a complete pinning of a Luttinger 
liquid with repulsive interactions, which results in a van- 
ishing conductance at zero temperature^*^ In addition 
to quantum wires, largely similar ideas apply to edge 
modes in a Hall bar geometry in the fractional quantum 
Hall regime, which are thought to behave as spatially 
separated chiral Luttinger liquids.'^ 

Evidence has recently emerged pointing towards the 
existence of the Luttinger liquid in metallic single-wall 
carbon nanotubes.^'^ The Luttinger liquid behavior was 
observed via the power-law temperature and bias- voltage 
dependence of the current through tunneling contacts at- 
tached to the nanotubes. Further technological advances 



have made possible the fabrication of low-resistance con- 
tacts between nanotubes and metallic leads (see, e.g., 
Refs. i8>9lilQlllll and references therein). These recent de- 
velopments have paved the way for systematic transport 
measurements in Luttinger liquids with impurities. 

In this paper, we study electron transport through 
a double barrier in a ID liquid. In the ID geometry, 
two impurities in effect create a quantum dot inside the 
system. Resonant tunneling through the two impuri- 
ties is a particularly attractive setup to investigate the 
correlated transport in an inhomogeneous Luttinger liq- 
uid. Due to the resonant behavior of the current, the 
interplay of Luttinger-liquid correlations and impurity- 
induced backscattering should be more easily accessible 
to transport measurements. Two striking experimental 
observations have been reported recently. In Ref. IT^ 
a resonant structure of the conductance of a semicon- 
ductor single-mode quantum wire was attributed to the 
formation (with reduction of electron density by chang- 
ing gate voltage) of a single disorder-induced quantum 
dot. In Ref. Il3l two barriers were created inside a car- 
bon nanotube in a controlled way with an atomic force 
microscope. In both cases, the amplitude of a conduc- 
tance peak Gp as a function of temperature T showed 
power-law behavior Gp oc T~'^ with the exponent 7 no- 
ticeably different from 7=1. The latter is the value of 7 
expected in the absence of interactions (for a review, see 
Refs. ITilTsI) provided T lies in the range F < T < A, 
where F is the width of a resonance in the transmission 
coefficient and A is the single-particle level spacing. The 
width of a conductance peak w followed a linear temper- 
ature dependence w oc T in both experiments. 

On the theoretical side, resonant tunneling in a Lut- 
tinger liquid was studied previously in a number of 
papers i2iiSiiLi2ii2i22, In particular, the width F cx 
was showniSiiS to shrink with decreasing temperature. 
The exponent ae depends on the strength of interaction 



2 



and describes tunneling into the end of a semi-infinite 
liquid. The dimensionless peak conductance (in units of 
e^//i) obeys Gp ^ T /T in the above range of T, which 
indeed leads to a smaller value of 7 = 1 — a^- The re- 
duced exponent 7 reported in Ref. [13 was positive (and 
different for different conductance peaks, in the range 
7 ~ 0.5 — 0.8), whereas in Ref. Il3 the reported value of 
7 ~ —0.7 was negative. More specifically, in Ref. 13. the 
conductance as a function of the gate voltage showed cer- 
tain traces of periodicity characteristic to the Coulomb 
blockade regime^^'i^'^ Surprisingly, both the amplitude 
Gp and the width w were reported to vanish with de- 
creasing T, in sharp contrast to the noninteracting case. 
While such behavior is known to be possible for very 
strong repulsive interaction^S*iS the required strength 
of interaction would then be much larger than expected 
and indeed reported (see Refs. Ial7l l2al23 and references 
therein) in carbon nanotubes. Roughly a doubling (or 
even a larger factor) of the expected.§*i22i2^ exponent Oe, 
which is Ue ~ 0.6 — 1.0, would be necessary to fit the ex- 
perimental data. Although the modelin g--'- ^ of carbon 
nanotubes as a four-channel Luttinger liquid has its share 
of complications, the observations^^ appear to present a 
puzzle. 

It was suggested in Ref. that a certain novel mech- 
anism of "correlated tunneling" dominates over the con- 
ventional sequential tunneling for T ^ F, leading to a 
doubling of the exponent a, which might explain the ex- 
periment. In the subsequent works^^iS^ the basic ideas 
relevant to the resonant tunneling have been questioned. 
In particular, the lowest-order contribution to the res- 
onance peak conductance Gp for A ^ T F has 
been argued^iiS^ to come from processes of second or- 
der in the end-tunneling density of states. While being 
a characteristic feature of the cotunneling regime far in 
the wings of the conductance peak, this, however, dis- 
agrees with the sequential tunneling picture inside the 
conductance peakjiSiiS This also poses a problem with 
the persistence^ of perfect transmission through symmet- 
ric barriers in the case of weak interaction. Moreover, the 
main suggestion24 is that taking a finite-range (but falling 
off fast beyond this range) interaction changes the con- 
ventional picture completely, as compared to the stan- 
dard treatment of a zero-range interaction in the Lut- 
tinger liquid model. According to Ref. [23. higher-order 
tunneling processes in combination with the effects of a 
non-zero range of interaction dominate the dependence 
of Gp on T even for T ^ F. A nearly perfect agreement 
with the experiment datai^ has been claimed (for a range 
of interaction far smaller than the distance between the 
barriers). However, no explicit formula for the conduc- 
tance has been given in Ref. while we see no ground 
for the low-energy long-distance physics to be essentially 
different if the radius of interaction is made finite. 

In another recent attempt^S to explain the 
experiments^ the observed power law was attributed to 
the contact resistance. Namely it was noted that the 
resistance of tunneling contacts to the leads, -Rc, and 



the resistance of the quantum dot add up (if one applies 
Kirchhoff's law), so that the anomalous T dependence 
in Ref. might be explained if Rc ^ Gp^, where Gp 
is understood as the resonance peak conductance of 
the dot. This is a legitimate proposition, although the 
measured contact resistance in Ref. ^3 was reported to 
be relatively lowiiSi^ 

It is thus desirable to examine the resonant tunneling 
in a Luttinger liquid in a broad range of temperature 
down to T = and for various parameters of the barri- 
ers. Our purpose in this paper is to analyze transport 
through a double barrier of arbitrary strength, strong or 
weak, symmetric or asymmetric, within a general frame- 
work of an analytical method applicable to all these sit- 
uations. There are a variety of techniques to construct 
the low-energy transport theoryi^i2L22, The method we 
develop here is valid for weak interaction and is based on 
the renormalization group (RG) approach of Refs.'2^'3l 
which was applied earlier in a variety of contexts. 
One of the appeals of this kind of theory is that it allows 
one to treat weak and strong scatterers on an equal foot- 
ing, which is technically significantly less straightforward 
in the bosonization methodjSi Thus it would be of inter- 
est to apply this method to a ID mesoscopic interacting 
system with many impurities (previously, the problem of 
transport in a disordered Luttinger liquid was studied by 
perturbative in disorder methods based on bosonic field 
theories in, e.g., Refs. I^ Ib^SsI) . A double barrier is the 
simplest "many-impurity" system which exhibits effects 
essential to the physics of disordered interacting ID liq- 
uids. 

Very recently, a generalization of the RG approach (ini- 
tially developed for a structureless barriesSSi^) has been 
proposed in Ref. for the model of a single impurity 
with an energy-dependent scattering matrix. This model 
is relevant to the resonant tunneling of weakly interacting 
electrons through a double barrier and accounts properly 
for the interaction processes within an energy band of 
width A around the Fermi level. In particular, a non- 
monotonic behavior of the conductance peak as a func- 
tion of T, caused by left/right asymmetry of scattering 
amplitudes, was investigated in Ref . |3^ within the single- 
resonance model. The fermionic RG™ has been shown to 
be in accord with the results^iSiiS for the resonant tun- 
neling in Luttinger liquids, obtained earlier by different 
methods. 

However, to describe microscopically the spatial struc- 
ture of a system of two or more impurities, a more sys- 
tematic analysis is needed. Specifically, one has to de- 
velop an approach that would include contributions to 
the renormalized scattering amplitudes due to interaction 
processes involving energy transfers larger than A. Also, 
the single-resonance model~ does not describe contribu- 
tions from tunneling through multiple levels inside the 
quantum dot (multi-level quantum dots represent a typ- 
ical experimental situation). Finally, it would be inter- 
esting to study transport through weak impurities (and 
through strongly asymmetric structures) for which the 
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bare transmission coefHcient exhibits no pronounced res- 
onant structure. AU this adds to our motivation to study 
the resonant tunneUng through a double barrier by gen- 
eralizing the RG approach of Refs. 

The RG approach enables us to investigate in detail 
the resonant transport of weakly interacting spinless elec- 
trons. Within the fermionic RG approach, we confirm 
earlier resultsM6_'i9 obtained within bosonic field theo- 
ries. We examine the conductance through a double bar- 
rier for arbitrary strength and an arbitrary shape of the 
barrier, not restricting ourselves to the model of a single 
resonant level. In particular, we demonstrate the exis- 
tence of narrow conductance peaks for two weak impu- 
rities, which is in sharp contrast to the noninteracting 
case. We do not find any trace of the correlated tun- 
neling mechanism proposed in Refs. ll3iii24.2{L We also 
clarify the relationship between the RG method^^i^ and 
Hartree-Fock (HF) approaches. 

The paper is organized as follows. First, in Sec. Ill Al 
we briefly outline the fermionic RG approach to trans- 
port through a single impurity. In Sec. Ill Bl we discuss 
a HF treatment of transmission through a single impu- 
rity and show its inadequacy to the problem. In Sec. lIIII 
we turn to a double barrier. We start with a perturba- 
tive expansion in Sec. lIIlXl and derive the RG equation 
for transport through the double barrier in Sec. IIII Bl 
We then analyze contributions to scattering amplitudes 
from different energy scales, compared to the level spac- 
ing inside the quantum dot and the resonance width, in 
Sees. IIII nillll El In Sec. IIII Fl we concentrate on the case 
of weak impurities. Finally, in Sec. II VI we calculate the 
amplitude and the shape of the resonance conductance 
peaks. 



II. SINGLE IMPURITY 

A. Renormalization group: Basic results 

We begin with a brief description of transport through 
a single structureless impurity in the spirit of the RG 
approach22i^. Without interaction, the impurity is char- 
acterized by a transmission coefficient to and reflection 
coefficients r^o and r^o i from the left and from the right 
respectively (we put the impurity at the center of coor- 
dinates, a; = 0). Suppose that the energy dependence 
of the bare scattering matrix can be neglected far from 
the boundaries of an energy band {—Do, Do) around the 
Fermi level. The energy scale Dq serves as the ultravi- 
olet cutoff of RG transformations and, physically, is of 
the order of vp/d (throughout the paper we put h = 1) 
or the Fermi energy ep, whichever is smaller.-- Here d 
is the radius of interaction and vp is the Fermi velocity. 
Deep inside the band {—Do,Do), we linearize the energy 
spectrum around the Fermi level. The differential RG 
equationg ^^i^° read 

dt/d£ =-atR, dTL.R/dC = arL.Rl , (1) 



where C = ln(Do/|e|), the energy e is measured from the 
Fermi level, the transmission probability T — and 
R = 1 — T. The boundary conditions at £ = set the 
scattering amplitudes at their noninteracting values to, 
r^o, and r^o- Throughout the paper we consider spinless 
electrons, for which the interaction constant is 

a^iVf- Vb)/2TTVF , (2) 

where Vf and Vb are the Fourier transforms of a pairwise 
interaction potential yielding forward {Vf) and backward 
(Vf,) scattering. The forward scattering does not lead 
to transitions between two branches of right- and left- 
movers, whereas the backscattering does. We assume 
that a > 0. 

Note that, for spinless electrons, the interaction- 
induced backward scattering and forward scattering re- 
late to each other as direct and exchange processes, so 
that the backscattering only appears in the combina- 
tion Vf — Vb and thus merely redefines parameters of 
the Luttinger model (formulatedi*^ in terms of forward- 
scattering amplitudes only). In particular, the backscat- 
tering does not lead to any RG flow for a. For spinful 
electrons this is valid only to one-loop order 

It is also worth mentioning that for a point interaction 
Vf — Vb, so that Q = 0, hence for spinless electrons one 
has to start with a finite-range interaction. However, the 
RG flow for the scattering matrix ^ occurs for |e| < 
vp/d and is governed solely by the constant a. It fol- 
lows that on low-energy scales one can effectively consider 
the interaction as local, Vcff{x — x') = 2TTavpS{x — x'), 
and formally deal exclusively with forward scattering. A 
non-zero range of interaction for kpd ^ 1 can mani- 
fest itself only in the boundary conditions to Eqs. at 
|e| ~ Do = vp/d and therefore does not affect the sin- 
gular behavior of the renormalized scattering matrix at 
e — > 0. We assume that the Coulomb interaction between 
electrons is screened by external charges (e.g., by metal- 
lic gates, in which case d is given by the distance to the 
gates) and that a resulting a ^ 1. For a treatment of the 
unscreened Coulomb interaction, see, e.g., Refs. I1I37I 

Integration of Eqs. (Q) givea^^'^° 



The phases of the scattering amplitudes are not affected 
by the renormalization. Equations 1^ are equivalent to 
a one-loop renormalization, so that Eq. |j2Jl is valid to 
first order in interaction 0(a) in the exponent of the 
power-law scaling. As follows from Eq. whatever the 
initial values of Tq, at a > they all flow to the fixed 
point of Eqs. Q at zero transmission,'^ T = at e = 0, see 
Fig-E In the limits of a weak impurity (both Rq -C 1 and 
R ^ 1) and a strong tunnehng barrier (Tq <C 1), Eq. Q 
coincides with the RG results obtained by bosonization^^ 
provided a ^ 1. Equation JSJ gives the transmission 
probability for electrons with energy e at temperature 
T — 0. For finite T, the renormalization stops at |e| ~ T. 
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FIG. 1: Schematic behavior of the transmission coefhcient 
T(e) for weak interaction and a single weak impurity with a 
bare coefficient To — 1 in three different regions, according to 
Eqs. @,l|3J: (f) a small logarithmic correction (5T; (II) power- 
law scaling of the correction 5T; (III) power-law vanishing of 
T. The Fermi level is at e = 0. 

Beyond the microscopic scale Dq; it is instructive to 
introduce two more energetic scales that characterize the 
renormaUzation of the transmission coefficient by weak 
interaction, and (Fig. 

ln(Z?p/i?o) = , DrlD^ = R^'" . (4) 

The energy Dp defines the scale on which a perturbation 
theory in interaction breaks down. If |e| Dp, the inter- 
action requires a non-perturbative treatment. The scale 
Dp does not depend on Rq and is much smaller than Dq 
for a <^ 1. The energy D^ defines the scale on which 
a perturbation theory in the impurity strength breaks 
down. If |e| < Z?^, a weak impurity with Rq ^ 1 yields 
strong reflection, R ~ 1. Provided Ro ^ 1, the scales 
Dp and Dr are parametrically different and, for any a, 
Dr <C Dp. We will see in Sec. IIIII that the scale Dr is 
of central importance in RG theory for a double-barrier 
structure. 

To summarize this section, there are three different 
types of behavior of T with e, as illustrated in Fig. 

I- -Dp ^ |e| ^ -Do, a small logarithmic correction 5T ~ 
aRo hidej/ZJo) is perturbative in both interaction 
and the impurity strength; 

II. Dr < |e| < Dp, power-law scaling of the correction 
(5T — Ro(|e|/I?o)~^" is perturbative in the impu- 
rity strength but non-perturbative in interaction; 

III. |e| < Dr, power-law vanishing of T ~ (|e|/Z?r)^" 
is non-perturbative in both the impurity strength 
and interaction. 

If the interaction is not weak (i.e., a is not small), the 
regime I shrinks to zero while the scaling exponents in 



the regimes II and III become —2a 2[(1 -I- 2a) — 1] 
and 2a 2ae = 2[(1 + 2aY^'^ - 1], respectively,^ 

To make a connection with experiment, it is worth not- 
ing that in the Luttinger liquid model (Vf, — 0) for N 
channels with the same Fermi velocity and no interchan- 
nel transitions, the constant a is related to the exponent 
tte which describes tunneling into the end of the liquid 
(see Sec. nj as follows: 

ae = A^"i[(l + 2iVa)i/2 - 1] . (5) 

E.g., iV = 1 for spinless and N = 2 for spin-degenerate 
electrons in a quantum wire with a single mode of trans- 
verse quantization,^ A'^ = 4 in the Luttinger-liquid model 
of the armchair carbon nanotubei^Si^^ Note that a^ — a 
for any N in the limit of weak interaction (a ~> 0). For a 
discussion of the tunneling into the end of a multi-channel 
Luttinger liquid with non-equivalent channels see Ref. .38t 

B. Hartree-Fock versus renormaUzation group 

We now turn to a conceptually important point that 
is highly relevant to our calculation of the resona nt tun- 
neling in Sec. IIIII The basic idea of Refs. I29l30l was to 
relate the RG transformations ^ to the scattering of 
an electron off Friedel oscillations created by the impu- 
rity screened by other electrons. ■^^ Indeed, to first or- 
der in a, when the interaction-induced correction to the 
density distribution may be neglected, the HF correction 
to the scattering amplitudes is logarithmically divergent 
oc ln(£>o/|e|)- The essence of this singularity is a slow de- 
cay of the Friedel oscillations, which behave in ID, in the 
absence of interaction, as sin(2fci?|a;| -I- 4'f)- Here 
kp is the Fermi wavevector, |a;| the distance to the im- 
purity, and 4'F a constant phase shift. This observation 
suggests that a HF treatment of the problem would yield 
a series of leading logarithms In the HF problem, 

it would be sufficient to solve iteratively in a a nonlinear 
Schroedinger equation [Hq -\- UhfH^} — e)'0 = 0, where 
Ho is the noninteracting Hamiltonian. The static self- 
consistent scattering potential is Ui -f Uhf, where Ui is 
the bare impurity potential and the nonlocal HF term 
reads 

UHFix,x')= - V{x - x')p{x,x') (6) 
-f 5{x — x') j dxiV{x — xi)p{xi, xi) . 

Here V{x) is a potential of pairwise interaction between 
particles, p{x, x') = J2q i^*qW)i^q{^) is a sum over all oc- 
cupied states, and p{xi,xi) in the Hartree term is the 
total electron density which includes the Friedel oscilla- 
tions. Self-consistency requires that both the wave func- 
tions il)q and the scattering potential Urf be corrected on 
every step of the HF iterative procedure. In effect, Ref.lsol 
(see also Ref. [s^) suggests that the HF procedure, being 
carried out in the "leading-log" approximation, would re- 
produce Eq. In particular, the HF calculation'^'^ of 
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second order in a gives a correction ~ a^£^ and, ac- 
cording to Ref. i3Ql this HF correction does coincide with 
the one obtained from an expansion of Eq. (PJ up to 
terms ~ 0{a^). Following this logic, the RG ^ is es- 
sentially an effective method of solving the HF problem 
self-consistently. 

In fact, however, while the HF theory is indeed a 
"logarithmic theory" in the sense that it can be ex- 
panded in a series of leading logarithms the series 
is very much different from the RG solution (PJ, even 
to one-loop order. To demonstrate this, let us repre- 
sent HF solutions for waves ^'±(2;) which are incident 
on the scatterer from the left and from the right 
(— ) as ip±{x) = a±+(x)e+(a;) -t- a±-(a;)e_(x), where 
e±(a:) — exp[±i(fci? -I- e/vF)x], aiiu(x) are amplitudes 
varying slowly on a scale of kp for kF\x\ 3> 1 and satis- 
fying the scattering boundary conditions a++(— cxj) — 
a__(oo) = 1, a_| (—00) = a ^(00) = 0. Separat- 
ing the slow and fast variables, the HF equations for 
ef/IeI 3> fc_F|a;| 3> 1 are written in the leading-log ap- 
proximation as 

da±^+/dCx = {a/2)T]a±- , 

da±^^/dC, = (a/2)?7*a±,+ , (7) 

where Cx = hi{DQ\x\/vp) and the function ri{x) describes 
the envelope n{x) = \ri{x)\ /2Tr\x\ of the Friedel oscilla- 
tions. The latter are given for kF\x\ 3> 1 by 

— Y.f *Re{(a;+a^_)e-2^(^-+^W''-} . (8) 

Solutions of Eqs. {Tj) for x > and x < are matched 
onto each other by means of the bare scattering matrix. 
Equivalently, Eqs. ((TJ can be cast in the form 

ain77/a£^ =a//2 , dl / dC^ ^ 2a\r]\'^ (9) 

with 77 = a++a!j + a i-'^^-- ^-n^d ^ ~ J2fj.u k^iyp- The 

function I{x) is a contribution of states close to the Fermi 
energy to a smooth part of electron density. Equations 
(ITI),© can be solved exactly. The scattering amplitudes 
and the shape of the Friedel oscillations are then found 
self-consistently. However, the rather intricate HF solu- 
tion appears to be a qualitatively wrong approximation 
in our problem of impurity scattering in an interacting 
ID liquid. The difference between the RG and HF solu- 
tions is already seen for a weak impurity with Rq ^ 1. 
While the RG gives for |e| 3> Dr a power-law behavior 
of the reflection amplitude 

rL,R = rw.RO e"^ , (10) 
the HF expansion yields a geometric progression 

1 , , 

rL,R = rLo,RO z 77 • (11) 



Equation Hll() can be obtained also by summing up lad- 
der diagrams in the particle-hole channel. Note a pole 
characteristic to HF approaches. 

We thus see that even for weak interaction the HF 
solution correctly describes the scattering by a screened 
impurity only for |e| S> Dp, where it gives a small per- 
turbative correction ~ aC. At the next order in a, the 
HF expansion does generate a term ~ tlo^hqo^ in r^^/j 
but, in contrast to Ref. l3Cl the numerical coefhcient in 
front of it is a factor of 2 larger than in Eq. which 
signifies a breakdown of the HF approach beyond the 
simplest Born approximation. 

It is worth noting that in Ref. another HF-type 
( "non-selfconsistent HF") scheme was employed to sub- 
stantiate the RG procedure. As mentioned above, single- 
electron wave functions within the self-consistent HF 
approach obey ■0? = 0° + GqUHF'ipq, where ipq sjnd 
the Green's function operator (5° describe noninteracting 
electrons. Let us write the correction to the transmission 
amplitude of second order in aC as Stl^e) = —CtoRo{oi£)^. 
The self-consistent HF approximation gives C = To-l-l/2, 
different from C — (3Tn — l)/2 following from the RG ex- 
pansion Q . In Ref. I32I only Uhf is renormalized while 
Ipq remain unperturbed, i.e., tpg ~* V'g in the right-hand 
side of the above HF equation. This yields C — 1, which 
accidentally coincides with the RG result for a weak im- 
purity. However, for arbitrary Tq this approach can be 
seen to fail as well, already at order 0{a^C^). 

It follows that the RG (Q) in fact cannot be obtained 
from the above HF schemes, i.e., the RG is not a method 
of solving the HF equations. There are other scattering 
processes, not captured by the HF approximation, that 
at higher orders in a almost compensate the HF result: 
a comparison of Eqs. H1U|I . (|11|I says that the coefficients 
c„ in the corresponding sums X)n '^n('^'^)" different 
in the RG and HF expansions by a factor of nl. 

One can see that the HF approach misses important 
scattering processes also by noticing that Friedel oscil- 
lations in a Luttinger liquid decay in the limit of large 
\x\ ^ vp/Dr as more slowly^S*^ (for the repul- 

sive interaction) than in a Fermi liquid. The Bragg re- 
flection by a potential created by the Friedel oscillations 
can then be easily seen to yield a power-law behavior of 
the transmission coefficient only if a is put equal to zero 
in the exponent of the Friedel oscillations. For any a > 
the Bragg reflection leads to an exponential decay of the 
wave functions, ln|0| cx — |a;|", instead of a power law 
required by self-consistency and by Eq. 

We can summarize the above observations by stating 
that, although it is certainly appealing to think of the 
RG Q as being associated with scattering off Friedel 
oscillations, the mechanism of the RG is more compli- 
cated. In fact, in a closely related calculation of various 
correlation functions of a clean Luttinger liquid without 
impurities the inapplicability of the HF theory was re- 
alized long ago>^2iiMi it is the interaction in a Cooper 
(particle-particle) channel that interferes with the HF in- 
teraction. As a result, instead of the summation of a 
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HF ladder one has to use a much more involved par- 
quet technique In the present problem, the par- 
quet summation is equivalent to the RG Q. A rigorous 
microscopic justification of the RG ^ can be done by 
using Ward identities in a diagrammatic approach and 
will be given elsewhere4& Our purpose in Sec. IIIII is to 
extend the RG ^ to the case of two impurities. 

III. DOUBLE BARRIER 

Consider now two potential barriers located a,t x — 
and X = xq and let the distance xq be much larger than 
the width of each of them. Clearly, the spatial structure 
itself yields an energy dependence of the total (describ- 
ing scattering on both impurities) amplitudes t{e) and 
^L,fi(e), even if such a dependence may be neglected for 
each of the impurities, as was assumed in the derivation 
of Eqs. Specifically, without interaction the energy 
A = TTVp/xo gives a period of oscillations in the total 
scattering amplitudes with changing Fermi level. If the 
impurities are strong, A is the level spacing inside a quan- 
tum dot formed by the barriers and the period of resonant 
tunneling oscillations. It follows that an RG description 
of scattering off a double barrier requires a generalization 
of the R C^^i^° to the case when the bare amplitudes to{e) 
and rio,i?,o(e) are energy dependent. 

A question, however, arises if it is at all possible to con- 
struct the RG theory for a compound scatterer in terms 
of only total S'-matrix, as in Eqs. Put another way, 
the question is if total scattering amplitudes generated 
by RG transformations are expressed in terms of them- 
selves only. As we will see below, the answer depends on 
the parameter A/Dq. We recall that Dq in our problem 
is the smallest of two energy scales given by and vp/d, 
where d is the radius of interaction. If A ^ Dq, so that 
there are many resonances within the band {—Dq,Dq), 
the RG transformations generate more terms than are 
encoded in the total S'-matrix. This is a generic case we 
are interested in. On the other hand, if A 3> -Do (e.g., 
when d ^ Xq), one would have from the very beginning a 
model of a single impurity with no spatial structure but 
with possibly energy dependent scattering amplitudes. 
In that case, it is sufficient to deal with total amplitudes 
only (see Sees. IIIIDIinEll . For the case of a single res- 
onance this is a model studied, e.g., in Ref. 47 (for an 
exactly solvable case of the Luttinger liquid parameter 
g — 1/2) and Ref. 13^ (for a weak interaction, 1 — g ^1). 

A. Perturbative expansion 

Building on our knowledge of the single impurity case, 
we start the derivation of a double barrier RG with a 
calculation of perturbative in a corrections to a time- 
ordered single-particle Green's function G'^^(x, x') — 
—i{T^^{x)ipj^{x')), where x = {x,t) and /x, ^ = ± label 
two branches of right (-1-) and left (— ) movers. The bare 



Green's function is characterized by bare scattering am- 
plitudes io(e), ''LO,flo(e)- To first order in a, transforming 
to the (x, e)-representation, which is most convenient in 
the present context, we have a correction to Gfi^{x{, xi; e) 
(summation over branch indices assumed): 

/oo 
dxGf^f^' {x,Xi;e) 
-OO 

X E^/j.'(a;)Gi.'^(a;f,a;;e) , (12) 

where 

'Sfj^,y{x) = iavp deG^u{x,x;e) . (13) 

"'--Do 

In the Luttinger liquid model, only forward scattering 
due to interaction is present: we explicitly assume this 
in Eqs. I|12() . H13|I by writing the self-energy that de- 
pends on two branch indices only [backward scattering 
can be straightforwardly incorporated for spinless elec- 
trons, see Eq. As explained in Sec. Ill Al we also as- 
sume in Eq. (|12|) that the interaction is effectively short- 
ranged and write the self-energy as a spatially local quan- 
tity. Indeed, the perturbative logarithmic correction to 
G ^^{xf,Xi]t) comes from energies |e| < vp/d, i.e., from 
spatial scales larger than d. We start with the case 
d <^ Xq (we will return to the simpler case d ^ xq m. 
Sec. IIIIDIi . Moreover, in the following, we neglect the 
forward scattering of electrons belonging to the same 
branch: to one-loop order, such processes can be seen 
to yield only a non-singular renormalization of the bare 
parameters. The forward-scattering interaction that gen- 
erates RG transformations similar to Eqs. is that of 
electrons from different branches 4^ We thus formulate 
an impurity problem with only non-diagonal couplings 

S+_(a;) and E y{x). 

The scattering amplitudes are related to G^j^^ixi, Xi] e) 
for e > (we measure e from the Fermi level upwards) as 

rL{e) - zi;FG+_(xf,Xi;e)e-"(^'+^-)/''-U,,,,^_oo , 
rnie) = zi;i.G_+(a;f , x;; e)e"(^'-^')/"- U^^.^^o , (14) 
t{e) = iVFG++{xi,Xi;e)e-''-^^'-''''>^^-\^,^oo,.,^-oo , 

and similarly for e < by changing ± — > =p in branch 
indices of G^,y. We count the phases of the reflection 
amplitudes from x — 0. Since the integral over x in 
Eq. H12|) involves integration over the interval of x inside 
the dot, < X < xo, corrections 6t{e), SrL^ii{e) cannot 
be expressed solely in terms of the bare amplitudes to{e) 
and rLo,flo(e)- A closed set of equations can be written 
by introducing amplitudes to stay inside the dot 

A^,_^(e) = ivpG^^.^ix, X- e)e2^'^'^ (15) 

(we need only nondiagonal amplitudes A^^_^), ampli- 
tudes to escape from the dot, to the left or to the right, 

d±(e) = zi;FG^,±(xf,x;e)e-'^("'-^")U,^±oo , (16) 
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and, similarly, amplitudes to get into the dot from outside 

6±(e) = iVFG±^^{x,xut)e'<^^^-^\.,,^^^ . (17) 

The amplitudes (|14|I - H17|I are constrained by unitarity 
and, moreover, &J^(e) = dy_^(e) by time reversal symme- 
try which we assume throughout the paper. In Eqs. H15|l - 
()17|l . X lies within the dot. Without interaction, the am- 
plitudes H15|) ~ (|17|l do not depend on x. This property is 
preserved in the leading-log approximation. 

The first-order corrections to t{e) and rL(e) read 

m = -^ ^{L+(e,e') (18) 
2 J-D„ e - e' 1^ 

+e{-e')t{e)[rni€)rl{t')x.-,+rL{e)rl{e')]] , 
SrLie) = -^ r ^{L_(e,6')+e(6')^L(e') (19) 

^ J -Da e ^ e 

+0(-e')[i^(6)r|,(e')x.-e'+ri(6)r2(6')]} , 

and similar equations can be written for corrections to 
rii{e) and the amplitudes defined in Eqs. p5|) - (|17(l . Here 
0(e) is the step function and 

X, = exp(2^ze/A) . (20) 

The terms L^(e,e'), which correspond to the integration 
over X in (|12ll inside the dot, are given by 

L^(e,e') = hl{e)A+Ae')d^_{e){xe-e'-l) 

+ 5:(6)A_+(e')<(e)(l-X.'-e) . (21) 

The amplitudes to stay inside the dot satisfy the relation 
Af^-^{€) = —A*_^^^{—e) and it is useful to decompose 
them as 

A^^_^(e) = e{e)B,,^,{e) + e{-e)C ,,^,{e) . (22) 

We thus arrive at a closed system of perturbative equa- 
tions written in terms of quantities describing the dot as 
a whole, without directly referring to parameters charac- 
terizing two barriers. 

B. Renormalization group for a double barrier 

Generically (singular exceptions are discussed below), 
the above first-order corrections diverge logarithmically 
as e ^ 0, which implies that higher orders of the per- 
turbation theory are important. Naively, one could try 
to treat Eqs. (|18|l . (|19() and other perturbative equations 
self-consistently, i.e., not as expressions for small cor- 
rections (Si(e), etc., but as equations to be solved for 
t(e) = to(e) + 5t{e)^ etc., where 5t{e) is not necessarily 
small. However, this would not correctly describe the 
case of strong impurities, either the case of a dot formed 
by two barriers or even that of a single strong barrier. 
To see this in the simpler case of a single impurity, take 



the limit xq — > 0. The terms L^(e,e') and factors Xe-e' 
then drop out in Eqs. p8|) . H19|l and the self-consistent 
equations acquire the form 

t[e) = io(e) (23) 

° Do 

Me) - ^^oW + f (24) 
^ J -Do e 

As a comparison of Eqs. ((TH|l . ((T^ and if^ . lf^ shows, 
the latter do not describe two impurities since they miss 
the terms Lp(e,e') and the factors Xe-e'- Moreover, 
Eqs. H23I) . H24|I do not describe even a single strong struc- 
tureless (with no dependence of the bare amplitudes on 
e) impurity. This can be checked straightforwardly, e.g., 
by putting all reflection coefficients in the integrand of 
Eq. (|23|l equal to unity: in this limit the integration over 
e' in Eq. H23(l gives a geometric progression of logarithms, 
instead of reproducing Eq. ^ for t{e). Also, plugging in 
Eq. |(2Jl into the integrand of Eq. (|24|l can be easily seen 
to give a logarithmic divergence at e — > 0, instead of a 
power law. The logarithmic divergencies are character- 
istic to HE approaches discussed in Sec. IIIBI Equations 
(|23|l . (|24|) . however, do describe correctly a weak struc- 
tureless impurity, but only for |e| S> Dr-. 

We now derive nonperturbative amplitudes for a dou- 
ble barrier using an appropriate RG scheme. To account 
for the e dependence of the bare amplitudes, the deriva- 
tion of the RG from the perturbative results ((TH|l . ((T^ 
necessitates introduction of two energies, e and D. The 
latter is a flow parameter in RG transformations, i.e., an 
ultraviolet cutoff rescaled after tracing over states with 
energies e' in the interval |e'| e {D,Do). The renormal- 
ization stops at D ^ max{|e|,T}. 

The essence of the RG procedurc^^'^° is a perturba- 
tive treatment of contributions to the renormalized am- 
plitudes at energy e from all states with energies e' in the 
interval |e'| g {D,AD), starting from D = Dq/A, such 
that A > 1 but alnA < l.'"^- The RG equations thus 
differ from both the HE equations and Eqs. lf^ . (|^ in 
that all the amplitudes depend on D and, moreover, the 
HF-type integration over projected states with energies 
e' only goes over the interval \e'\ E {D,AD) instead of 
{0,Do). 

In effect, each step of the RG transformations accounts 
for the scattering off the Friedel oscillations in a finite 
spatial region, |a;|, jx — xq] S {vp / AD,vf/ D)- Moreover, 
the Friedel oscillations are only partly modified, through 
the (already performed) renormalization of the reflection 
amplitudes at energies larger than D. At the same time, 
the scattering matrix at energies smaller than D is taken 
at its bare value. This should be contrasted with the 
HF approach, where the scattering amplitudes are de- 
termined by interaction processes on all energy scales on 
every step of the HF iterations. 
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The system of one-loop RG equations for a double bar- 
rier reads 



where 



dt(e,D) 
dCn 



dCn 



h>{e,D)[L+{e,e';D) 

e{-e')t{e,D)[rn{e,D)r*j,{e',D)xe-e' 
rUe,D)rl{e',D)]} , (25) 

L{e,D) [L-{e,e';D) + 9{e')rLie,D) 



[{e,D)rl{e\D}]} , 



(26) 



and similar equations for other amplitudes. Here Cd — 
\ii{Dq/D) and we introduced D dependent amplitudes 
t{e,D), etc., so that integration of Eqs. over D 

acts on all the amplitudes [which should be contrasted 
with Eqs. (|23|l . (124(1 . where the integration acts on only 
one amplitude in the products of three] . All amplitudes 
in Lfj^{e, e') H21|l are now also functions of D. The integral 
operator I^' (e, D) is defined as 



Ie'{e,D) 



a 




r] 


de' 


21nA 


J D 


J-AD 


e-e' 



7{---} : (27) 



where A ^ 1 is restricted by the condition a In A ^ 1. 
The theory is renormalizable if the action of the operator 
/e' (e, D) to leading order is independent of A, which is 
the case for the present problem within the leading-log 
approximation. Needless to say, in the limit a;o ~* [i-e., 
L±{e,e';D) and Xe-e' 1] Eqs. 123,® describe 
a single impurity. 

At finite temperature T, one should substitute the 
Fermi distribution function np{e) for the step func- 
tions in Eqs. (|25|) . H26I) and also in Eq. (|22|l according 
to 0{±e) npiT^)- The factor (e - e')~^ in Eq. ^ ef- 
fectively stops the renormalization at D ~ |e|, while the 
factors np{±e') do so at I? ~ T, otherwise the renormal- 
ization can be carried out down to D = 0. The infrared 
cutoff at _D T establishes a characteristic spatial scale 
of Lt — vp/T. Due to the thermal smearing, the Friedel 
oscillations decay exponentially on a scale of Lt [which 
can be seen from Eq. ||SJ) if one plugs in np{e) into the 
integrand and extends the integration to e > 0] . The ab- 
sence of renormalization at I? ^ T is closely related to 
the thermal cutoff of interaction-induced corrections^^ to 
the conductance in higher dimensions. 

The RG equations (|2Sl,(|211) should be solved with 
proper boundary conditions at _D = Dq, where the renor- 
malization starts from the bare values of the amplitudes 
((T^ - l(T7l) . For the double-barrier problem, the boundary 
conditions for the transmission and reflection amplitudes 
are: 



^o(e) 



tlt2 

S{e) 



fLo(e) = ^"1 + 



r2ti 
S{e) 



(28) 



Sie) 



(29) 



and riio{e) = -r2g(e)to(e)/i(5(e) by unitarity. The coef- 
ficients ri^2i'ri 2) E^re the noninteracting reflection ampli- 
tudes from the left (right) and ti^2 are the transmission 
amplitudes of each of the two barriers, respectively. Sim- 
ilarly for other amplitudes: 



d+ie) 
dZ{e) 



r[ h 
d+(e) ^ to(e) 

f 2 Xe t2 

toU)r2 



(30) 



Xe , C+_(e) = - 



We count the phases of ^ from x — Q and the phases 
of ri_2 are defined for an impurity sitting at x = 0. 



C. Separate renormalization of two 
impurities: D ^ A 

We are now in a position to solve the system of RG 
equations ((231,(1201) by integrating out all states with 
energies |e'| > max{|e|,r}. We begin with the case 
Dq ^ A, which is a typical case unless interaction is 
very long ranged. We proceed in two steps. Let us first 
integrate over D 3> A. This can be done for arbitrary e. 
Specifically, if |e| > A, this will already solve the problem 
by providing us with fully rcnormalized amplitudes. In 
the more interesting case of |e| < A, we will only sum up 
contributions to the rcnormalized amplitudes from states 
with |e'| > A and, as a second step, will have to proceed 
with renormalization for < A. 

Since the renormalization for D 3> A involves many 
resonant levels, the amplitudes contain slowly varying 
parts and parts oscillating rapidly with changing e' on a 
scale of A. Integration over e' in Eq. H27|l allows us to 
separate the slow and fast variables: as a result, the de- 
pendence of the amplitudes on D will be slow on the scale 
of A. However, even after that the RG equations look 
rather cumbersome. To construct the solution to these 
equations, note that an important parameter D / D^-^^^ is 
available, where 



-^r-min = minjUj-j^, 13^2} 



(31) 



and Dr^ 2 are defined for each of two barriers by Eq. Q). 
If both barriers are initially (i.e., dX D = Dq) strong 
(1^1,2! ^ 1), then this parameter is small for all D < Dq. 
However, if one or both of the barriers arc initially weak, 
there is a range of 13 G {Dr^-^, Dq) where at least one 
barrier still remains weak. 

It is useful first to examine some general properties 
of integrals over e' that appear in the course of renor- 
malization. Let us return to the perturbative expansion 
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(|18|l . (|19|) . We see that the averaging over e' involves two 
types of integrals 



Sie') 



To 



(32) 



where Sit) is given by Eq. (|29l) and Xi^2 are related by 
X2 = [Ii — C)/r'ir2- The integrals (pl^ are evaluated in 
different ways depending on whether at least one of the 
barriers is weak (so that \r'ir2 \ <C 1) or both barriers are 
strong (|r']^r2| — 1). In the former case the integrand of 
Ii is only slightly modulated, so that one can expand the 
factor S'~"'^(e') and average over harmonics x"(£')- Then 
only zero harmonics give rise to singular (logarithmic) 
corrections and in the leading-log approximation we have 



Ii^C , I2 







(33) 



In the opposite case of strong barriers, sharp resonances 
appear that are described by a Breit-Wigner formula for 
(e') and yield |Zi| ~\l2\: 



Ii = C/2 , I2 



-C/2r[r2 



(34) 



The situation repeats itself in the RG equations H25|I . H26|I . 
The difference in the factor of 1/2 between the values 
of Xi in Eqs. H33() . H34(I implies that the renormalization 
should be carried out differently (see Fig.|21) in the regions 
D > i?r„i„ and A < £) < -Dr„i„. 

For D ^ i'r^in, similarly to Eq. (|33|l . after the averag- 
ing over e' only zero harmonics contribute to the renor- 
malization. The solution at D ^ -Dr-mi,, is then simple: 
it has the form of Eqs. (|28|l . (|30|) with the reflection and 
transmission amplitudes of each of two barriers renormal- 
ized separately (Fig.|2^), according to the RG CJ,© for 
a single impurity: 



dti,2{D) 
dri,2{D) 



= ^ati,2{D)\n.2{D)f , (35) 
= an,2{D)\h^2{D)\^ ■ (36) 



One can check straightforwardly that Eqs. I|28|l . (|30|l (de- 
scribing the Fabry-Perot resonance) with the replacement 



H,2 



ri.2iA.2) 



\T^.2V + {D/D^Y-\h^2?Y'^ 
^l,2(?-i^2) 

\ri,2? + {D/D^Y-\h^2?V'^ 



(37) 
(38) 



solve Eqs. l|^ .lf ^ averaged over harmonics for D ^ 
max{Z?rj^;^, A}. 

On the other hand, if I?r,„i„ A, there is an interval of 
D e (A, Dr^^J) in which each of the impurities is strongly 
reflecting (Fig. I^la), so that the averaged equations are 
again simplified by summing over resonance poles, simi- 
larly to Eq. (|34|l . We get an independent renormalization 
of ti^2{D) according to 



h.2{D)/tl,2iDr^J = (^/A,„.„)"^' 



and the scaling exponent is now half that for D Dr^^^ ■ 
We thus have two solutions given by Eqs. (|S7|l . and 
Eq. (|39|l . respectively, that match onto each other at I? ^ 
Z^r^in . Due to the slow power-law dependence, for a ^ 1 
the matching is exact. 




(39) 



FIG. 2: Sketch of different stages of the RG procedure for a 
double barrier for i3 > A > T, Eqs. The Friedel 

oscillations shown schematically around each of the barriers 
yield a renormalization of the scattering amplitudes and are 
renormalized themselves as well. At energy e, the renormal- 
ization comes from spatial scales smaller than = «F/|e| 
(these regions are marked by the solid lines), (a) Separate 
renormalization of two weak barriers. Since there is no en- 
ergy quantization between the barriers, the Friedel oscillations 
both inside and outside the dot contribute to the renormaliza- 
tion and the barriers "do not talk to each other" . The scaling 
exponent a is the same as for a single barrier. Sec. Ill AI (b) 
Separate renormalization of two strong barriers. This figure 
describes either initially strong barriers or those that are ini- 
tially weak but both become strongly reflecting due to the 
renormalization. The inner part of the dot with a discrete 
energy spectrum does not contribute to the renormalization. 
The renormalization of each of the barriers is governed by the 
exponent a/2, twice as small as for weak barriers. 

We conclude that the key difference between the renor- 
malization for D larger and smaller than Dr^^^ is that 
for D ^ ^r„in the transmission amplitudes for each 
barrier are renormalized with the exponent a, whereas 
for D -C Dr^in with the exponent a/2. In both limit- 
ing cases, two barriers are renormalized separately for 
Z? S> A (Fig. 121). It is worth stressing that generally the 
independent renormalizations of two barriers cannot be 
derived from RG equations written in terms of only t(e) 
and rL_ij(e), i.e., the terms L^{e,e'; D) are of crucial im- 
portance in the derivation of Eqs. (|35|l " (|38|l . However, 
the renormalization of resonant tunneling amplitudes for 
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energies near resonances allows for another formulation 
which involves t{e) and jj(e) only, we will return to 
this issue in Sec. IIII Dl 

If e is close to one of resonant energies, Eqs. (|25ll . (|26() 
can be further simplified by expanding Xe near the res- 
onance: the renormalized amplitudes for strong barri- 
ers take then the form of Breit-Wigner amplitudes with 
D dependent widths Ti^2{D) = (A/27r)|ii,2(I?)l^ oc D°'. 
Specifically, for initially strong barriers: 



(40) 



where |ti.2p ^ 1 ai'e the bare transmission probabil- 
ities at D — Do and resonant peaks are sharp [i.e., 
ri,2(-D) A] for all D < Dq. If at least one barrier 
is initially weak, the resonant structure develops only at 
D ^ Dr^i^- Provided one barrier is initially weak (as- 
sume this is the right barrier and 0^-2 3> A), whereas the 
other is strong, then 



(41) 
(42) 



If A ^ S ^ Do, i.e., both barriers are initially 
weak, then 



boundary conditions in the double-barrier case should 
be written at D '--^ A, instead oi D ^ Do- We obtain for 
kl,k'l,^«A: 



dt{e,D) 
drL{e,D) 



-- t{e,D)[rR{e,D) 7%iD) 
+rL{e,D)TliD)] , (44) 

-rl{e,D)TliD)] , (45) 



and similarly for ri^(e, D), where the bar over the reflec- 
tion amplitudes denotes the averaging (|27|l over e'. 




(43) 



and r2{D) is given again by Eq. (|42|l . 

To summarize this section, substituting D — s- |e| in 
Eqs. (|37|I . (|38(I . H39|I and using the Fabry-Perot equations 
()28|) gives the fully renormalized scattering amplitudes 
for |e| > A if |e| > T. Also, if T > A, substituting D -> 
T solves the problem for arbitrary e. However, when both 
\e\,T <t; A, we should proceed with the renormalization 
in the range D ^ A. 



FIG. 3; Illustration of the renormalization of a double barrier 
for D < A, Eqs. The Friedel oscillat ions are shown 

schematically within the range = vf/\^\ outside the dot, 
where they contribute to the renormalization of the scatter- 
ing amplitudes at energy e (provided |e| S> T). The double 
barrier for small D may be considered as a single barrier with 
an energy dependent scattering matrix, which describes the 
resonance (sketched by the dotted line) on a single level. The 
scattering matrix inside the renormalized resonance is dis- 
cussed in Sec. IIII El 



D. Single resonance: D <^ A 

Let us now consider Eqs. for |e|, |e'| < A. 

In this limit, the terms H21I) containing the amplitudes 
A^,_^(e') to stay inside the dot become irrelevant in the 
RG sense: the phase factors Xe-e' in Eq. (PT|l can be ex- 
panded about e, e' = 0, which leads to the cancellation of 
the singular factor (e — e')~^ in Eq. (|27|l . As a result, the 
terms L^t(e,e';I?) do not contribute to the renormaliza- 
tion at I? <C A. The factors Xe-e' should also be omitted 
in the terms of Eqs. (|25|l . (|26|) that are proportional to 
r|j(e'). Thus we are led to a coupled set of RG equa- 
tions that describe also a single impurity with energy 
dependent scattering amplitudes: the spatial structure 
of the double barrier system is of no importance for the 
renormalization at Z) ^ A (see Fig. However, the 



We integrate now Eqs. (|44|) . 145|) assuming that each of 
two barriers is characterized by Dr^ j 3> A. This condi- 
tion means that either the barriers are strong initially at 
D — Do or get strong in the course of renormalization 
(|35|l . (|36|) before D equals A. We will analyze the case of 
both or one of Dr^ ^ being smaller than A in Sec. IIII Fl 

Thus, if 2 3> A, we have narrow peaks of resonant 
transmission in a neighborhood of the Fermi level, even 
if the bare reflection coefficients are small. Consider first 
the case of a resonance energy eg lying exactly on the 
Fermi level, eo = 0. Let r(D) be a renormalized width 
of the resonance peak at the Fermi energy (to be found 
below). If L» > V[D), then |rT:^(£))| ~ 1, which allows 
for a significant simplification of Eqs. (|44|l . (|45() . It is 
convenient to introduce phase-shifted amplitudes — 

—iu} ~ — -(-27ri(eF+eo)/A r , —i(u}t +iOt 1 

r^e ^^^,rii = rfie '2 , t = te HVti+vta^^ 
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where ipr-^ is the phase of ri, etc., in obvious notation. 
Then we get, by putting the averaged ampUtudes far from 
the resonance fL,R{D) = 1: 

= tie, D) [f^ie, D) + fn{e. D) ] , (46) 

" i [ ''L.iile: D) -*(£,£))], (47) 



with the foUowing solutions 



rL{e,D) 



4(I3)-m1(D)]1/^ 

u+{D) + 2ie 
u_(D) +2ie 
u+{D) + 2ie ' 
-u^{D) + 2ie 
u+{D) + 2ie ' 



(48) 
(49) 
(50) 



Ti>T2 




where 



FIG. 4: Transmission coefficient T(e) for temperatures Ds <^ 
T <^ A. Two curves correspond to two temperatures: Ti 
(solid) larger than T2 (dashed). While the peak height is 
not renormalized by interactions, the resonance gets narrower 
with decreasing T, Eq. H52^ . 



Ui{D)^T±iA){D/Ar, (51) 

and r±(A) = ri(A) ±r2(A) should be found by match- 
ing onto Eqs. I|40|I ~ H43() . The width of a resonant tunnel- 
ing peak r{D) is thus given by u+{D). 

Note that the only condition we have assumed in the 
above derivation is Z? ^ r(Z?) with D — max{|e|,T}, 
otherwise e in Eqs. H48|I - H5U|) may be arbitrary. Thus, 
Eqs. l|^ - (|5n|l give the shape of the e dependence of fully 
renormalized amplitudes for the case of temperature T 3> 
r(T) (with T substituted for D). In particular, Eq. ^ 
says that the width of the resonance behaves as T": 

r(r) = r+(A)(r/A)°. (52) 

As follows from Eq. (|48|l . while the resonance becomes 
sharper with decreasing T, the peak value of the trans- 
mission amplitude is not renormalized (see Fig. ^ , since 
the D dependent factors cancel in Eq. H48|) at e = 0. The 
absence of renormalization stems from the vanishing of 
the sum fi(e, D) + fnie, D) in Eq. at e = 0, which 
can be seen from Eqs. 

We recognize Eqs. (jI5|) - (|5(H) as Breit-Wigner solutions 
that take into account renormalization at Z? < A. On the 
other hand, we have already obtained Breit-Wigner for- 
mulas in Sec. IIII CI where the renormalization has been 
carried out for £) A, for e close to a resonance energy. 
In particular, the results of Sec. IIII 01 applv for |e| ^ A 
if Co = 0. The matching of the two solutions at 13 ~ A 
implies that Eqs. (|44|l . (|45|) are in fact valid in a broader 
range of D, namely for D -C Dr^i„, provided only that 
one averages rii{e',D) over e' together with the phase 
factor Xe'- It follows that in the case of strong barri- 
ers close to resonances the RG equations can be cast in 
the form (|^ . (|T7l) containing t{e, D) and rL,R{e, D) only. 
Note also that for Dq < A the boundary conditions to 
Eqs. l|I^ . (|T7|l are fixed at _D = _Do, which leads to the 
change A ^ Dq in Eq. (|FTT) . 



At this point, one might be concerned about a possible 
contribution to t(e = 0, D) from other resonances. In- 
deed, in the derivation of Eq. 148|) , which gives no renor- 
malization of t{e = 0,13), we approximated |r(e',£')| 
by unity for large \£'\,D ^ r(_D). Corrections coming 
from other resonances arc clearly small in the parameter 
^i,2{D)/ A but one should check if they might contribute 
to the renormalization of t{e — 0, D). Relaxing the above 
approximation by allowing for resonant "percolation" of 
electrons through the barriers at _D > A does give a per- 
turbative correction to the RG (|46|l : 

— = -at{e,D) [rL{e,D) -rR{e,D)] , 

(53) 

which, in contrast to Eq. H46|l . does not vanish at e = 
(unless the double barrier is symmetric: the correction 
is then always zero). However, Eq. H53I) tells us that the 
correction is irrelevant since U-{D) itself scales to zero 
as D" . We thus conclude that the "single-peak approx- 
imation" of Eqs. (|46|l ■ (|47|) correctly describes the renor- 
malization of the resonant amplitudes for all D ^{D). 

E. Inside a peak: D ^ Ds 

Now that we have integrated out all D ^ u+{D), let 
us continue with the renormalization for D inside a res- 
onant tunneling peak. The point at which D and m+ (D) 
become equal to each other yields a new characteristic 
scale Ds- 

Ds - r+iA){Ds/Ar = r+(A)[r+(A)/A]"', (54) 

where r+(A) is obtained from Eqs. 140|I - H43() depending 
on the ratio of 2 ^^'^ ^o- To leading order in a the 
exponent a' = q;/(1 — a) — > a. As will be seen below. 
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the significance of Ds is that the width of the tunnehng 
resonance saturates with decreasing D on the scale of Dg- 
For D <s: Ds, the RG equations (03,(03 can be sim- 
pUfied since the scattering amplitudes now depend on 
a single variable, which is D — max{|e|,r}. The av- 
eraged reflection amplitudes rLj^{D) coincide then with 
fL.R{D), and the RG equations can be written in pre- 
cisely the same form as for a single impurity: 



dt{D) 



-at{D)R{D) , 



drL,RiD) 



(55) 

with matching conditions aX, Cd — ^"i^iDa/ Dg). The dif- 
ference between the single structureless impurity and the 
resonance peak is that in the latter case the ultraviolet 
cutoff is Dg. The fact that the scattering amplitudes in- 
side the peak are described by Eqs. H55|l was used also in 

Ref . m ____ 

In the preceding Sees. IIII GIIII'DI asymmetry of the 
double barrier, i.e., a possible difference between ti and 
t2 was seen to determine the amplitude and the width 
of a resonance peak but otherwise did not lead to any 
qualitatively different consequences, as compared to the 
symmetric case. However, for small D <^ Dg the renor- 
malization of the scattering amplitudes is essentially dif- 
ferent depending on whether the barriers are identical to 
each other or not. 

Consider first the symmetric case. Assuming, as in 
Sec. IIII 01 that the resonance energy eg = 0, we get the 
Breit-Wigner formula with a width given by Eq. (|51|l : 



m = 



u+{D) 



u+{D) + 2ie 



(56) 



which turns out to be valid down to Z? = 0. A remarkable 
consequence of Eq. 1)56(1 is that the resonance in the sym- 
metric case is perfect, T = 1 at e = 0. While this is trivial 
for noninteracting electrons, weak interaction in the Lut- 
tinger liquid is seen to preserve the perfect transmission, 
in agreement with the result obtained by a bosonic RG.^ 
So long as inelastic scattering is not taken into account, 
the perfect transmission at e = is not affected by finite 
temperature T, either, as can be seen from Eq. (|56|l if one 
puts D — T . However, the width of the resonance does 
depend on T. At T = 0, the width is finite and given by 
Dg (see Fig.|3Ji), which follows from Eq. H56|) for D = |e|. 
For T:^ Dg, the width obeys Eq. lf3^ . 

The shape of the perfect resonance depends on the 
parameter T/Dg. If T ^ Dg, the reflection probability 
as a function of |e| behaves near e = first as for 
|e| < T and then as |e|2(i-") for T < |e| < Dg. For 
larger energies, the transmission probability falls off with 
increasing |e| as 



T(e) = (i^,/2|e|)2(i-«) . 



(57) 



This lineshape should be contrasted with the Lorentzian 
which describes the transmission peak for T 3> Dg up to 
|e| ~ T, at which point a crossover to Eq. H57|) occurs. 



Let us now turn to the asymmetric case. Inspecting 
Eqs. (|55|l . we see that a new characteristic scale Z3_ 
emerges: 



' ' r+(A) ) 



l/a 



(58) 



which coincides with Dg for strongly asymmetric barriers 
but vanishes for symmetric ones. For T > D- we get the 
same results for T(e) as in the symmetric case, only with 
an overall factor of 



A = 



r2+(A)-r2_(A) _ 4ri(A)r2(A) 



n(A) 



[ri(A) + r2(A)] 



(59) 



However, for T <C -D- a new feature in the behavior of 
T(e) shows up, namely a power-law falloff with decreasing 
|e|. The function T(D), as obtained from Eqs. ((55|l . for 
D <^ Dg reads 



T{D) = 



X{D/Dg 



l-X[l-{D/DgY-- 
One sees that T(e) behaves as (Fig. [S^) 

T{e)^\{\e\/D^f" 



(60) 



(61) 



in the interval T <^ |e| <^ _D_ and saturates at smaller 
energies: T(e = 0) = X{T/D-f". Thus, in the hmit 
T <^ D^ , the resonant transmission probability as a func- 
tion of e exhibits a double-peak structure, see Fig.lSJ^. If, 
however, the barriers are only slightly asymmetric, the 
gap near e = develops in a range of e which is much nar- 
rower than the width of the resonance peak. Specifically, 
T(e) first grows up wit increasing |e| for T ^ |e| <C D^, 
then there is a plateau with an energy independent trans- 
mission for D^ ^ |e| -C I?s, and T(e) starts to fall off as 
e| is further increased. 

In the above, we analyzed the behavior of T(e) for the 
resonance energy eo — 0, i.e., when it coincides with the 
Fermi level. Let us now examine the case cq 0. Again, 
let the barriers first be symmetric. Then, in Eq. H56I) , the 
resonant denominator changes to u+{D) + 2i{e — eq) and 
D is, as before, max{|e|,T}. The innocent looking shift 
e — > e — Co leads at T = to dramatic consequences for 
transmission at the Fermi energy, e = 0. Namely, T(e) 
is now seen to vanish at e = and zero T, whatever cq 
unless it is exactly zero. A new characteristic scale Di 
becomes relevant at eg ^ 0: it is defined by u^{Di) = 
2|eo|, which is rewritten as 



D,=Dg{2\eo\/D, 



(62) 



The significance of the energy Di is that the width of 
the gap in the dependence of T(e) around e = at T = 
is given by Di for |eo| < Dg. Note that Di <C |eo| for 
|eo| < Dg. 

The shape of the resonant peak as a function of e 
changes in an essential way for |eo| < Dg. Specifically, 
if T 3> the changes are weak; however, for T Di 
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gap width = D ^ 

(b) 



FIG. 5: Schematic summary of transmission peak structure 
for zero temperature and strong barriers: (a) eo = 0, a single 
peak of T(e) for symmetric barriers with a width Ds trans- 
forms into a double peak with a reduced height and a gap 
width D- < Ds for asymmetric barriers, Eqs. H58|I . H6()|I : (b) 
< |eo| < £^s, a single peak for symmetric barriers, centered 
at £ = £0, shows a dip at the Fermi energy, e = 0, with a 
gap width Di < |eo|, Eq. 1621 . For £o / and asymmetric 
barriers, (a) describes the case of sufficiently small £o, namely 
D- S> Di, whereas (b) with a properly rescaled height de- 
scribes the opposite case. The only effect of finite T < Ds is 
to fill the gaps at |£| < T. 



a range of e arises, T ^ |e| ^ Di, within which T(e) 
behaves as (Fig. 



T(e) - (|6|M) 



2a 



(63) 



The power-law fallofF (|63|l saturates close to the Fermi 
level at T(e = 0) = {T/Dif°'. 

We thus see that the width of the resonance in the 
transmission through a symmetric barrier exactly at the 
Fermi energy T(e = 0) as a function of eg vanishes as 
T ^ 0. On the other hand, the width of the resonance in 
the transmission at eg = a-s a function of e is finite even 
at T = and is given by Dg. This peculiar feature is in 
sharp contrast to the resonant tunneling of noninteract- 
ing electrons, for which the two widths are the same. 

For asymmetric barriers, T(e) does not change sub- 
stantially with increasing |eo| as long as Di -C and 
is given by the formulas for symmetric barriers with an 
overall reduction of the transmission probability by a fac- 
tor of A (|59l) otherwise (see Fig.l^J. 



rescale [according to Eqs. H40|l - H43|) ] parameters in oth- 
erwise general formulas for the resonant tunneling. Let 
us now examine the resonant transmission in the case 
of at least one barrier being initially so weak that the 
renormalization does not make it strong at _D ^ A. 

We begin with the case of both barriers characterized 
by 2 ^ The total reflection coefficient rL{e,D) 
as obtained from Eqs. ()25fl . H26|l in the limit [r^l ^ 1 is 
given by 

rL{e,D) = {n-r2Xe){Do/Dr (64) 

with D = max{|e|,T}. It is worth noting that, due to 
the oscillating factor Xe (I20|l - it is not possible to derive 
Eq. (|64(l from Eqs. (|23|l . (|24() even in the simplest case of 
a weak double barrier. 

Suppose first that the barrier is symmetric and eo = 0. 
Then Eq. simplifies to 

R(e, L») = 2[ 1 - cos(27re/A) ] {Dr/Df°'. (65) 

One sees that reflection is enhanced by interaction, but 
the reflection coefficient is always small, R <C 1 for any e, 
provided Dr ■€i A. No sharp features in the e dependence 
of the scattering amplitudes emerge around the Fermi 
energy (see Fig.EJi). 



T(e) 



T(8) 



sym 



asym 




All 
(a) 



FIG. 6: Schematic behavior of the transmission coefficient for 
zero temperature and weak barriers, Eq. 17H : (a) for eq — 0, 
the modulation of T(£) for both symmetric (solid) and asym- 
metric (dashed) barriers is strongly increased. However, for 
symmetric barriers T(£) remains close to unity for all e, while 
there appears a dip of width 5- < Dr, Eqs. 16711 . 16811 . around 
£ = for asymmetric barriers; (b) for £o 7^ and symmet- 
ric barriers, the bare transmission (dash-dotted) is strongly 
renormalized (solid) due to interactions and exhibits a gap 
of width Si < Dr, Eq. 17UII . around the Fermi level. For 
£0^0 and asymmetric barriers, (a) describes the case of suf- 
ficiently small £0, namely 5- 3> (5i, whereas (b) with a prop- 
erly rescaled height describes the opposite case. Similarly to 
Fig.lHl finite T fills the gaps at |£| < T. 



F. Weak barriers 

In Sees, nil (^III El we have assumed that -Dr„i„ ^ A, 
which means that even if the impurities are initially 
(at D = Do) weak, they get strong in the process of 
renormalization before D becomes equal to A. Under 
this assumption, we have sharp resonant peaks close to 
the Fermi level and the bare scattering amplitudes only 



It is now instructive to introduce a weak asymmetry 
R_ = IR2 - Ri|, such that R_ < R ~ Ri^2- Given that 
the asymmetry is weak, it can manifests itself only at 
small energies. Expanding Eq. H64|) about e = 0, we get 
for |e| < A: 



R{e,D) = 
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As can be seen from Eq. H66() . asymmetry sets two new 
characteristic scales of energy: (R_/R)A and a smaller 
scale 



= Dr{R-/2R) 



l/a 



(67) 



Provided that temperature T <C (R_/R)A, the reflec- 
tion coefficient starts to grow with approaching the Fermi 
level at |e| ^ (R_/R)A. The enhancement of reflection is 
cut off by temperature before R becomes of order unity 
if T is not too low, specifically if (5_ <C T. However, if 
T <C (5-, then reflection gets strong at |e| ^ To de- 
scribe the scattering probabilities at |e| < (5-, one should 
solve Eqs. ()55|l . derived in the same way it was done in 
Sec. IIII El now with matching onto the perturbative (in 
Ri^2) solution (|64|l anywhere in the region S- ^ |e| ^ A. 
For D < (R_/R)A the solution reads: 



T{D) = 1/[1 + {S^/D) 



2q1 



(68) 



Thus, the weak double barrier remains slightly reflecting 
after the rcnormalization provided that T ^ S-. How- 
ever, for both T,\e\ <^ the transmission probability is 
small: within the range T ^ |e| <C (5-, T(e) behaves as 
(Fig.Et) 



T(6) = {\e\/S^) 



2a 



(69) 



and saturates for smaller |e| at T(e = 0) = {T/S^)^°'. 
Comparing Eqs. H69|l . (|61|l with each other, we see that 
the energy 5- is a counterpart of D_ for the case of a 
weak barrier. 

Generalizing to ep 7^ 0, we have for |eo| ^ A a shift 
e ^ e — eo in Eq. (|66|) . while D ~ max{|e|,r}. A new 
energy scale <^ eo appears, at which the reflection 
coefficient becomes of order unity in the symmetric case: 



(5i - I?, (2^|eo|/A)i/", 



(70) 



analogous to D\ in Eq. H62I) for tunneling barriers. The 
energy 8\ gives the width of the gap in the transmission 
probability at the Fermi level at T = 0. For T < |e| <C 
(5i, we get a power-law vanishing of T(e) = (|e|/5i)^" with 
decreasing |e| (see Fig. I^Jj) and a saturation for smaller 



at (r/(5] 



\2a 



A general expression for the scattering 



probabilities, valid for arbitrary D^^ ^ ^ A and |e|, D <C 
A, can be obtained from Eqs. H55(l : 



R(e,g) 
T(e,C) 



(r} 



/2 



R, 



1/2 



(71) 



(RiR2)^/Mf ) 



27r 



D 



Equation H71(l reproduces Eqs. (|66|l - l|7()(l in the corre- 
sponding limits. 

We conclude that if the barriers are symmetric but eo 
is nonzero, or if the barriers are asymmetric, the trans- 
mission probability vanishes (Fig. ^ at the Fermi level 
in the limit T ^ 0. We will see in Sec. IIVI that these 



features lead to the emergence of a sharp peak in the 
low-temperature conductance as a function of eg even for 
two weak impurities, provided only that they are slightly 
asymmetric. 

Finally, when two strongly asymmetric barriers are lo- 
cated nearby, so that at -D A one barrier is strongly 
reflecting whereas the other is still weak, the effect of the 
latter on the transmission probability remains small for 
any D. Let us take the example ^ A and 13^2 ^■ 
Then we get for Z?,.^ > Z5 > A 



T{e,D) = {D/Dr.f' [l + 2(D,J£))"cos0] 



(72) 



where 6 = 27r(e — eo)/A. One sees that the presence of 
the weak impurity only leads to a weak modulation with 
changing e. For _D < A, the independent rcnormalization 
of the weaker impurity is suppressed by reflection from 
the strong barrier and T(Z5) behaves as (D/Z^n)^" down 
to L» = 0. 



IV. CONDUCTANCE PEAK 

The solution to the problem of transmission through a 
double barrier given in the preceding sections allows us 
to examine the linear conductance of the system G{eo, T) 
as a function of temperature T and the energy distance 
between the Fermi level and a resonance level eo- Recall 
that we have studied the elastic transmission of inter- 
acting electrons, i.e., the energy e of an incident elec- 
tron before and after the transmission is the same, while 
in the process of scattering off the barrier e is not con- 
served due to interaction with other electrons both in 
the particle-hole and Cooper channels. At flnite T, there 
are also inelastic processes, characterized by the inelastic 
scattering length ijn. Neglecting the inelastic scattering 
is legitimate if ^ Lt ~ vp/T^ which is satisfied in 
the present problem for weak interaction a <^ 1. On 
the other hand, it is worthwhile to note that the very 
formulation of the scattering problem in the interacting 
case even for elastic scattering is a delicate issue if one 
keeps scaling exponents of higher order than linear in 
a.^^ Also, interaction-induced current-vertex corrections 
in the Kubo formula in the presence of impurities and 
interaction are governed by exponents of higher order in 
a. Here, we avoid these complications by treating the 
scattering problem within the one-loop approximation, 
i.e., keeping only first order terms in the exponents. Un- 
der these conditions, one can use the Landauer-Biittiker 
formalism relating the conductance and the transmission 
probability. The conductance G'(eo,T) in units of /h 
reads 



G(eo,T)= deT{e){~dnFlde) 



(73) 



We are interested in the low-temperature regime with 
T <^ A, otherwise we intend to keep T arbitrary, i.e., T 
may be as small as zero. 
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Below we analyze various limiting cases. As we have 
seen in Sec. Illll the strength of interaction, the strength 
of the barriers, and the degree of asymmetry set a num- 
ber of characteristic energy scales which yield a variety of 
different regimes in the temperature and energy depen- 
dence of the transmission probability T(e,T). By means 
of Eq. H73|l . these regimes manifest themselves in the be- 
havior of the conductance G'(eo,T), which is a directly 
measurable quantity. 

Before turning to the calculation of G(eo,T), let us 
briefly discuss the restrictions and possible implementa- 
tions of our model. Even though we are dealing with the 
case of weak interaction, a <^1, the product aln(I?o/A) 
may be large in small quantum dots, so that the renor- 
malization of scattering amplitudes on scales D ^ A is 
necessary, as has been done in Sec. IIII CI Experimen- 
tally, as discussed in Sec. Q the parameter a is typi- 
cally not small in carbon nanotubes [in Refs. 'gIJ the 
value of a extracted according to Eq. |(SJ from the mea- 
sured ae ~ 0.6 — 1.0 is a ~ 1.3 — 3.0)]. On the other 
hand, in single-mode semiconductor quantum wires the 
strength of interaction is typically smaller (in Ref. 
ae 0.2 - 0.5 corresponds [Eq. ®] to a ^ 0.3 - 0.8) 
and can be more easily made small through screening 
by nearby metallic gates. In weakly interacting wires, 
our results would be applicable directly. At the same 
time, our theory captures much of the essential physics 
of strongly [with a defined in Eq. |2l of order unity] cor- 
related wires too, as follows from the comparison with 
the known results obtained by different methods Ai&iS 

As for the impurity strength, both transport and tun- 
neling experiments on quantum wires with imperfections 
have been so far focused on the case of initially strong in- 
homogeneities (impurities, artificially created tunneling 
barriers, or non-adiabatic contacts). Sections IIII CIIII'EI 
describe this situation in detail. Moreover, as shown in 
Sees. IIII CIIIIFI weak inhomogeneities, which are poten- 
tially realizable in a controllable way both in semiconduc- 
tor quantum wires and carbon nanotubes, are of special 
interest, since the Coulomb interaction transforms two 
weak impurities into a quantum dot with a pronounced 
resonance structure. Finally, a relatively strong asymme- 
try of the quantum dot appears to be inevitably present 
in the experimental setups of Refs. 1121.131 . 



A. Strong barriers 

Consider first the case of strong barriers (more pre- 
cisely, the bare transmission through the barriers may 
be high, T1.2 — 1, but we assume that the barriers 
get strong before the RG flow parameter D equals the 
single-particle energy spacing inside the dot. A), i.e., 
Dr^^^ = Do(min{Ri,R2})^/2" > A. Then we have a 
sharp peak of the transmission probability centered at 
e — eo whose width is max{Z?s, r(T)} <C A, where Dg 
and r(T) are defined in Eqs. (|5^ . (|^ . In other words, 
the width of the peak in T(e,r) is r(T) = DsiT/D,)" 



for T ^ Dg, whereas for smaller T <^ Ds the width is of 
order Dg and does not depend on T . 



1. T ^ Ds, Sequential tunneling 

For T ^ Ds, the shape of the conductance peak is 
given by (see Fig. 

where the peak value of the conductance 

Gp = 7rAr(r)/8r , (75) 

with A defined in Eq. and C = (max{|eo|, T}/T)". 

The width of the conductance peak w is of order T, as 
for noninteracting electrons; however, the power-law be- 
havior of Gp{T) is seen to be modified by in teract ion, in 
accordance with the results derived in Refs. ll6ll9i Note 
that the scaling of Gp oc T""^ is governed by the single- 
particle density of states pe (T) for tunneling into the end 
of a Luttinger liquid, namely Gp oc p^(T)/Tm^ 

We recognize Eqs. (I74|) . (|75|l as the conventional se- 
quential tunneling formulas, but with a T depen- 
dent resonance width T{T). Far in the wings of 
the resonance the exponential falloff H74(l . G(eo,T) ~ 
Ar^^r(|eo|) exp(— |eo|/r), crosses over onto the cotun- 
neling (determined by the processes of fourth order in 
tunnehng amplitudes) power law G(eo, T) = Ar^(r)/4e§, 
as usual. The crossover between the sequential tunneling 
and cotunneling regimes occurs at |eo| — Tln[T/r(T)]. 
In fact, this formula is valid^^ for an arbitrary strength of 
interaction with r(T) oc Pe{T) oc T°'' , where ae (equal to 
a for a weak interaction) is the end-tunneling exponent 
(O. It is worthwhile to note that for strong enough in- 
teraction (namely for ae > 1) the sequential mechanism 
of tunneling is effective for the resonance peak for all T, 
down to T = 0>i^ Moreover, for ae > 1 the crossover 
to the cotunneling regime shifts towards larger |eo| with 
increasing strength of interaction. 

As demonstrated in Ref. IT^ Eqs. l|71|l . (|75|) can be ob- 
tained from a classical kinetic cquation'^"'^ for occupation 
numbers characterizing the state of the quantum dot. We 
get the same results from the fermionic RG equations. It 
is worth mentioning that, although Eqs. (|74|l . (|75|) follow 
from a classical kinetic equation^ ^ which involves diago- 
nal elements of the density matrix only, the transmission 
in our derivation is fully coherent. The fact that in the 
high-temperature limit it can also be described in terms 
of the classical kinetic equation merely means that quan- 
tum corrections to the kinetic equation may be neglected 
for small T(T)/T [whereas the quantum suppression of 
the tunneling density of states Pe{T) may be treated in 
this equation in a phenomenological way] . 
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FIG. 7: Strong barriers: Qualitative dependence of the con- 
ductance G (in units of e^/h) on eo, the energy difference 
between the Fermi energy and the resonance level, for dif- 
ferent temperatures T and weak interaction, (a) Resonance 
contour for symmetric barriers. The peak conductance Gp 
grows [Eq. (I75ll 1 with decreasing T and saturates [Eq. 1761 1 
at Gp — 1 as T — > 0. The width of the resonance shrinks 
to zero, according to Eqs. 17411 . 1771 1. (b) Resonance contour 
for asymmetric barriers. The dependence of Gp on T is non- 
monotonic: Gp grows [Eq. 17511 . curves T4 and T3] with de- 
creasing T for large T and goes down [Eq. 1791 . curves T2 
and Ti] for T < D-. The width of the resonance decreases 
[Eq. with lowering T and saturates [Eq. (731] as T -» 0. 



2. T <^ Ds, Symmetric barriers 

Let us now turn to low temperatures T -C -Ds, where 
processes of all orders in the tunneling ampUtudes are 
important. Consider first the symmetric case (Fig. EK)- 
The main contribution to the integral over e in Eq. H73|) 
comes from |e| ~ T <^ Dg, so that the shape of the 
conductance peak is a Lorentzian: 



G(eo,r) 



r2(r) 



(76) 



We see that the height of the peak Gp — 1 and the width 



w = Ds{T/DsT 



(77) 



exhibits a power-law temperature dependence with an 
exponent depending on the strength of interaction. The 
vanishing of w as T — > should be contrasted with the 
behavior of the peak in T(e, T) (Fig. (SJ, whose width is 
Ds for low T. In the limit T — > 0, the conductance peak 
becomes infinitely narrow but the resonance at en = 
persists down to T = 0, in accordance with Ref. We 
thus confirm the persistence'^ of the perfect resonance at 
T = by means of the fermionic RG. For finite T <C Ds, 
there is a small correction 1 — [T / Dgf'^^~"'\ which 

comes from the non-perfect transmission for finite e at 
£0 = after the thermal averaging ifT^ . 



3. T <C Ds, Asymmetric barriers 

We recall that the renormalization in the asymmetric 
case is governed by the scale D_ [defined in Eq. (|58|l ] 



which describes the degree of asymmetry. The double- 
peak structure of the transmission coefficient as a func- 
tion of eo for T ^ D_ translates into a complete vanish- 
ing of the conductance peak at T — > 0. Specifically, for 
T :^ the conductance G(eo,T) is given by Eq. l(7H|) 
for symmetric barriers with an overall factor of A. How- 
ever, for T <^ D_ the transmission at the Fermi level 
falls off with decreasing T (Fig. 13)) and so does the con- 
ductance peak (see Fig. [3d): 



G(eo,T) 
which gives 

Go — A 



A(T/D, 



\2a 



(D_/Ds 



(2eo/Ds)2 



(78) 



T 

dI 



r+(A) ' 

where r±(A) are defined below Eq. H51|) . Thus, at small 
T ^ D_ , the height of the conductance peak goes down 
as T decreases, whereas the width of the peak does not 
depend on T any longer. This kind of behavior of the 
resonance peak was predicted in Ref. 5 Gp{T) scales as 
p^iT), as for a single impurity. The role of asymmetry 
was emphasized and expressions similar to Eqs. (17911 w ere 
obtained in Ref. 1^3. However, the authors of Ref. l33l do 
not distinguish between the scales D_ and Ds (the latter 
is denoted e there). 



B. Weak barriers 

Consider now the case of weak barriers, i.e., D^^ ^ <^ 
A. Naively one could think that scattering on weak bar- 
riers cannot possibly yield a sharp peak of conductance. 
Indeed, the transmission probability as a function of e 
(Fig. 13 does not have any peak at e = eo, in contrast to 
the case of resonant tunneling. At high T ^ D^, G(eo, T) 
is a weakly oscillating (with a period A) function of eo. 
The only difference with the non-interacting case is an 
enhanced amplitude of the oscillations. 

Let us show that in fact the interaction-induced van- 
ishing of T(e) at the Fermi energy e = for T = does 
lead to a narrow Lorentzian peak of G(eo, T) (see Fig.lHJ, 
provided that T is low enough and the barriers are not 
too asymmetric. 



1. Symmetric barriers 

Integration l|73|) of the transmission probability H71|) 
for symmetric barriers yields 

-1 



G(eo,r) = 



1 



2^eoy ,^D,^2" 
T 



(80) 



which indeed describes a Lorentzian peak (Fig. ISji) with 
the height Gp = 1 and the width 



A / T 



TT VD, 



(81) 
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It follows that the peak is narrow, w <C A, provided that 
T ^ Dr- In the limit T — > 0, the width of the peak is 
infinitesimally small. Similarly to the case of resonant 
tunneling, the resonance at ep = remains perfect in the 
presence of interaction at T = 0. At finite T, there is a 
small correction 



goes down for w ^ T <C A, where the correction behaves 
similarly to the case of symmetric barriers, 1 — Gp cx 
y2(i-Q) rpj^g conductance peak is narrow provided the 
asymmetry is weak, R_ <C R. If the asymmetry is strong, 
R_ ~ Ri + R2, the peak is completely destroyed. 



l-Gp^iT/Af{Dr/T) 



2a 



(82) 



Equation H82|l describes also the high-T behavior of Gp 
in the case of slightly asymmetric barriers (see below). 




A/2 



FIG. 8: Weak barriers: Qualitative dependence of the con- 
ductance G (in units of e^/h) on the resonance energy eo 
for different temperatures T and weak interaction, (a) Sym- 
metric barriers. The value of G(eo — 0) grows [Eq. 1821 1 
with decreasing T and saturates [Eq. 1811 1 at Gp = 1 as 
T — > 0. The width of the resonance shrinks to zero, accord- 
ing to Eqs. 1801 . 1811 . The resonance peak becomes visible 
at T ^ Dr (curves T2 and Ti). (b) Asymmetric barriers. 
The dependence of G(eo = 0) on T is non-monotonic: it 
grows [Eq. 1821 1 with decreasing T for large T and goes down 
[Eq. H84|l ] for T < 6-. The width of the resonance decreases 
[Eq. 1811 . similarly to the symmetric case] with lowering T 
and saturates [Eq. I|84|l . curves T2 and Ti] as T — > 0. If asym- 
metry is strong, G(eo) is only slightly modulated for all T, 
i.e., exhibits no peak. 



2. Asymmetric barriers 

Introducing a weak asymmetry R_ = |Ri — R2I ^ R ~ 
R-1,2, we get for T <^ S-, where 6- is defined in Eq. (|^ : 



G(eo,r) 



R^{T/Dr 



R^/4-f R2(27reo/A)2 
The height and the width of the peak are 



(83) 



2a 



Thus, the asymmetry leads (Fig. ^p) to vanishing Gp at 
T ^ and the width is seen to saturate with decreasing 
T, similarly to Eq. H79() . It is worth noting that the 
dependence of Gp on T is non-monotonic: Gp oc T^" 
grows with increasing T for T ^ (S_ , continues to grow in 
the range <^ T ^ w according to 1 — Gp oc T~^°', but 



V. CONCLUSIONS 

In conclusion, we have thoroughly studied transport 
of weakly interacting spinless electrons through a dou- 
ble barrier. We have described a rich variety of differ- 
ent regimes depending on the strength of the barrier, its 
shape, and temperature. We have developed a fermionic 
RG approach to the double barrier problem, which has 
enabled us to treat on an equal footing both the reso- 
nant tunneling and resonant transmission through weak 
impurities. In the latter case, we have demonstrated how 
the interaction-induced renormalization in effect creates 
a quantum dot with tunneling barriers with a pronounced 
resonance peak structure. Moreover, we have shown that 
even very weak impurities, for which the renormalized 
transmission coefhcient does not exhibit any peak, may 
give a sharp peak in the conductance as a function of 
gate voltage, provided that the double barrier is only 
slightly asymmetric. In contrast, the resonant structure 
is shown to be completely destroyed for a strongly asym- 
metric barrier. All the regimes we have studied may be 
characterized by three different types of behavior of the 
conductance peak height Gp and the peak width w on 
temperature T: 

(i) for high temperature T 3> r(r), 
Gp oc T"^^ and w ccT; 

(ii) for lower T, depending on the shape of the barrier 
(whether it is symmetric or asymmetric), either 

Gp does not depend on T and w oc T" or 

(iii) Gp oc T^" and w is constant. 



One can see that none of the regimes (i-iii) supports 
Gp oc T^""^ and w oc T, as proposed in Ref.lla. Further 
experiments would be useful to resolve the puzzle. In- 
cluding spin and generalizing to the case of several chan- 
nels (possibly with different Fermi wavevectors) within 
the framework of the present approach warrant further 
study. 
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